On the Newtonian origin of the spin motive force in ferromagnetic atomic wires 
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We demonstrate numerically the existence of a spin-motive force acting on spin- carriers when 
moving in a time and space dependent internal field. This is the case for electrons in a one- 
dimensional wire with a precessing domain wall. The effect can be explained solely by adiabatic 
dynamics and is shown to exist for both classical and quantum systems. 

PACS numbers: 75.47.-m, 72.25.Rb, 75.75.+a 



I. INTRODUCTION 



Recently Barnes and Maekawa 1 have proposed a gen- 
eralization of Faraday's law to account for a non- 
conservative force of spin origin. This arises in sys- 
tems with time-dependent order parameters as a result 
of Berry phase (BP) accumulation 2 . As an example they 
consider a domain wall (DW), formed in a finite ferro- 
magnetic wire and precessing about a static co-axial ex- 
ternal magnetic field. In the adiabatic approximation, 
where electron spins remain aligned with the local mag- 
netization, a constant potential shift Acf) is generated be- 
tween the two ends of the wire. This is directly propor- 
tional to the angular frequency of precession of the wall 



A(j) = — uj , 
e 



(i) 



and within the Stoner model it exactly cancels the Zee- 
man potential. Such a potential, described as a spin- 
motive force (SMF), has been recognised previously in 
the context of the Aharonov-Casher 3,4 and Stern's 5 ef- 
fects. These are all manifestations of BP related phe- 
nomena, where holonomies arise as a result of a parallel 
transport of some kind 6 . The latter does not need to 
be a quantum effect, another example being the classical 
Foucault pendulum. 

Here we demonstrate computationally the result of 
Eq. (1) through time-dependent quantum-classical simu- 
lations of an atomic wire incorporating a precessing DW. 
We also present an analytical classical argument for the 
driving mechanism of the SMF in this system. Our ap- 
proach has the benefit of being "Berry-phase- free" in the 
sense that it does not need to call for a Berry phase argu- 
ment to explain the SMF and demonstrates the Newto- 
nian nature of the conversion of the magnetic response of 
electronic spins into an electrostatic voltage drop. This 
is further illustrated with classical dynamical simulations 
for a system of classical magnetic dipoles in a rotating 
magnetic field mimicing the DW. In addition we show 
that if one abandons the Stoner model and accounts for 
a non-spin component of the magnetic moments forming 
the DW, the cancellation between the SMF and the Zee- 
man potential is incomplete, leaving behind a non-zero 
net SMF, which can be experimentally measured. 



II. MODEL 

We consider a one-dimensional magnetic atomic wire 
in a magnetic field and describe the conduction electrons 
by an s-d 7 tight-binding Hamiltonian 



where (cf ) is the creation (annihilation) operator for 
an electron with spin ±1/2 (a = 1, 2) on the atomic site 
i and cr is the vector of Pauli matrices. The first term in 
Eq. (2) is the spin-independent tight-binding (TB) part, 
while the second describes the spin interaction with the 
effective local field 7 , 
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where Eq is the onsite energy (Eo = OeV for all sites), 
n = e 2 /47T£o = 14.4 eV A, \ is the hopping parameter, g e 
is the electron ^-factor and B is the external magnetic 
field. The second term in the brackets in Eq. (3) is a 
mean-field repulsive electrostatic potential 8 with an on- 
site strength U and a Coulombic decay at large intersite 



distances Rij. Aqi = qi — qp is the excess number of 

electrons on site z, q^ being the equilibrium one. 

In Eq. (4) is the effective local angular momen- 
tum at site z, associated with the magnetic moment of 
that atom, normalized by h (and thus dimensionless). 
are treated as classical variables, nonetheless exchange- 
coupled with strength J > to the conduction electrons 
according to a classical Hamiltonian 

Us = - J] Si • [J(a)i + 5 sMbB] - J z (Si • z) 2 • (5) 

i i 

Here #s is the ^-factor of Si which could be of mixed 
spin and orbital origin, z (|z| = 1) is the unit vector 
and J z the anisotropy constant along the easy z-axis. 
(s)i = Tr [pucr] is the expectation value of the electron 
spin at site z, where p is the density matrix and the trace 
is over the spin coordinates. 

The corresponding quantum and classical Liouville 
equations of motion for the two subsystems are 



dp 

dt 



p,n e 



^ = {S Z ,H S }, (6) 



2 



(a) 

m 



1 1 


Is 













FIG. 1: (Color online) Different prospect views of the DW 
formed by the local spin {S^} in the middle of a single-atom 
chain. An external magnetic field, applied along the wire, 
induces a clockwise rotation of the DW about the z-axis. 



where {•, •} represents the classical Poisson bracket. 
These have been solved numerically with a quantum- 
classical dynamical simulation following the evolution of 
a preset DW structure in a metallic atomic wire in mag- 
netic field. Initially the set of classical spins {S^ ) is 
prepared in a DW arrangement (see Fig. 1) and relaxed 
self-consistently in the electronic environment. At time 
t = the external magnetic field B = Bz is switched on 
along the wire and a new initial electronic state is self- 
consistently determined for {Si}(°). The system is then 
propagated according to Eqs. (6). 



III. NUMERICAL SIMULATIONS 

The electrostatic potentials Vl(r), developing away 
from the DW on the left (right )-hand-side of the chain, 
are computed as the spatial (over two identical sets L(R) 
of Ny atoms at each wire end) and temporal (over the 
evolution time T) averages of the onsite potential, i.e. 

t n 



VL(R)(T) = 1/(TN V ) 



KAq n (t) 



where 



N is the total number of atoms. We investigate the sta- 
tionary voltage drop AV ca i c = lim T — oo [Vl(T) - Vr(T)] 
that builds up across the system. In the limit of local 
charge neutrality (U —> oo) this should compensate for 
any arising spatial energy disturbances in the system. We 
anticipate two such contributions 



AV = - g e fi B B/e 



(7) 



where the first term is due to the proposed non- 
conservative SMF from Eq. (1), while the second is due 
to the Zeeman split. In order to extract the effect of 
the SMF itself in the first set of simulations we have set 
g e = 0, so that only the first term in Eq. (7) remains. 

The parameters used for the simulations are x — 
— 1 eV, qf^ = 1.75 e/atom for every z, U = 7 eV, 
J = 1 eV, J z = 0.5 meV, g s = 2, a = R^ i+1 = 2.5 A. 
TV = 400 atoms so that the chain is much longer than 
the typical width of the relaxed DW (about 10 atomic 
spacings). The values of x, J and are chosen such as 
to produce a halfmetallic system with a completely filled 
spin-up band, which lies about 0.5 V below the Fermi 
level. 

The spin-DW in these simulations undergoes a steady- 
rate clockwise rotation with an angular frequency oo 
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FIG. 2: (Color online) Time evolution of some dynamical vari- 
ables at B = 100 T and for g e = 0: (a) S^ and S y local spin 
components at the DW center zo , showing the clockwise rota- 
tion of the DW about the z-axis. The angular frequency cj of 
the DW precession is extracted by fitting S x (T) to cos (ooT) ; 
(b) longitudinal displacement of the DW center zo; (c) aver- 
aged potentials Vl, Vr and Altaic (see text). 
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FIG. 3: (Color online) Calculated SMF as a function of the 
DW precession and dependence of the slope over the Coulomb 
parameter U for g e = and gs = 2. (a) The calculated 
stationary AV depends linearly on u with a slope h ca \ c ~ 
0.92/1 for realistic values of the parameters J and U; (b) h ca \ c 
tends to saturate at the exact value of h with increasing U. 



about the direction of the field [see Fig. 2(a)] and ex- 
hibits nearly rigid oscillations about a center zq, slightly 
displaced to the left [see Fig. 2(b)]. The steady ro- 
tation generates a SMF, with gentle oscillations that 
correlate with those of the DW (since the projection 
of the total spin in the system on the direction of the 
field is conserved) and which has an asymptotic time- 
averaged value AV^aic [see Fig. 2(c)]. The dependence 
A Kaic (a;), obtained by sweeping the external field be- 
tween 20 T and 500 T, is linear [see Fig. 3(a)] with a slope 
^caic = 0.606 eVfs ~ 0.92 h. The deviation of h ca \ c from 
the exact h [from Eq. (1)] is studied with respect to the 
two main assumptions in our model: (i) the adiabaticity, 
which is governed by the strength of the exchange cou- 
pling J, and (ii) the local charge neutrality, which allows 
us to identify \AV\ with the SMF and is exact only for 
U — > oo. The first criterion is found to be well satisfied 
for J = 1 eV. Increasing J ten times results in less than 
1% improvement in h ca \ c . The ratio h ca \ c /h, however, is 
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found to be sensitive to U and it asymptotically tends to 
1 as U is increased [see Fig. 3(b)]. This result confirms 
the validity of Eq. (1) and indeed demonstrates that a 
SMF originates from the precession of the DW. 

In reality, however, the effect of the applied magnetic 
field on the electrons cannot be switched off. We, there- 
fore, return to Eq. (7) and rewrite it in the form 



at/ h 9e H 
AV = —co ous 

e gs e 



9e \ ft 
1 -I -co. 

9 s 



(8) 



Here cos — gsl^BB/ft is the Larmor frequency of the lo- 
cal spins. The actual angular frequency of precession of 
the DW co differs slightly from cos due to the exchange 
interaction with the conduction electrons. In order to 
account for this effect, we have introduced an effective 
#g such that co = g^jiBB/ft. We have verified Eq. (8) 
numerically by varying the value of gs (see Fig. 4). The 
effective value g$ is determined by the calculated preces- 
sion frequency of the wall (Fig. 4c). Finally we have 
again obtained a value of ft ca \ c ~ 0.92 ft, identical to the 
previous finding in the case g e — for this choice of ex- 
change parameter and charging strength. 
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FIG. 4: (Color online) Computational demonstration of Eq. 
(8) for g e — 2 and a set of values gs = 0.5 -r- 3. Panel (a) 
shows the linear dependence of the stationary potential drop 
Altaic on the angular precession frequency co; (b) is used to 
determine the effective g-factors g$ and they are compared to 
the input values gs in (c). Note that g$ = gs for gs = g e = 2. 
Panel (d) demonstrates the validity of Eq. (8). 

Apparently, the voltage drop across the system fully 
disappears when gs = #g = g e -> as derived in reference 
[1] for the Stoner model. However, in s-d systems, where 
#g has a partially orbital origin, this is not the case and 
the SMF manifests itself as a measurable quantity. This 
could be used to determine the effective ^-factor of the 
localized spins. In particular if the DW precession is 
blocked, the measured drop would be just equal to the 
Zeeman split, i.e. a measurement could determine if the 
wall is precessing or not. In the remaining part of the pa- 
per we examine the mechanism for the SMF by a classical 
analogy. 



IV. CLASSICAL PORTRAIT 

Instead of quantum electrons as in the first part of the 
paper we now consider non-interacting classical particles 
with an intrinsic angular momentum s (|s| = s = ft/ 2) 
and with the electron mass m e . These are trapped in 
a one-dimensional box, where an inhomogeneous time- 
dependent field 

b(z, t)=b (cos(/t) sin(6y , sin(/t) sin(0*), cos(0*)) (9) 

is present and 6 Z = 6(z) is chosen such as to mimic a 
continuous DW-like structure 9 , rotating rigidly with an 
angular frequency /. Here / is analogous to co from the 
quantum simulation, though / > corresponds to an 
anticlockwise rotation about the longitudinal axis. The 
classical Hamiltonian of the spin-particles in the field b 
is analogous to that of the quantum electrons interacting 
with local spins {S^}, 



Hdass = 7> 7 s ' b (^, t) , 



(10) 



where 7 is the coupling strength (analogous to J) and p 
is the canonical momentum of the particles. Then Hamil- 
ton's equations of motion 10 are 

m e z = 7s • X7 z h(z, i), s = 7s x b(z, t) (11) 

where V z = (d/dz) Sj t- 

We consider the limit of large 76, in which the dynam- 
ics of the spin-particle becomes adiabatic, in the sense 
that s remains closely aligned with b and its precession 
about b is by far the fastest motion in the system. How- 
ever, for s to follow b(z, £), there must always be some 
residual misalignment 11 between the two. This is nec- 
essary in order to generate those torques, which, when 
averaged over the quick precession of s, enable s to keep 
up with h(z,t). This small misalignment, marked by the 
angle cp in Fig. 5, is also the origin of the effective New- 
tonian force on the spin-particle that manifests itself as 
SMF. 

Fig. 5 depicts s and b at some instance of the par- 
ticle's migration. Differentiating the relation between 
the angles at the bottom vertex of the tetrahedron 
cos a = cos up cos — sin up sin 9 cos ft under the condition 
Scpi = s sin up sin ft = const (which corresponds to keeping 
s and t fixed), in the adiabatic limit up — > 0, we obtain 



v z ip = - cos (ft) v z e. 



(12) 



From Eq. (11) the longitudinal force F z and the torque 
T = I T I = 75^6 are related by 

F z = - 7 |s| |b| sm(up) V z up = -TV z up = T 2 V Z Q , (13) 

where Eq. (12) has been applied and the full torque T de- 
composed into two orthogonal torques with magnitudes 
Ti = Tsin^ and T 2 = T cos ft (see the inset of Fig. 5). 

In the adiabatic regime (s^ <C 5), we average the two 
components of the torque over the fast precession of s 
about b. These averaged torques T\ and T 2 must be driv- 
ing the two separate motions of the spin as the particle 
crosses the region of the rotating DW-like field, namely, 



FIG. 5: (Color online) A snapshot of spin-particle's passage 
through the DW-like region of h(z,t). 
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a rotation in the b-z plane enabling s to keep up with the 
spatial variation of b, and another rotation in a plane 
perpendicular to the z-axis which makes s follow the an- 
ticlockwise precession of b, and thus 

7T«siV*0, -|s x f| = -sf sm9. (14) 

Applying Eq. (13), the averaged linear force upon 
the spin-particle in the rotating magnetic field is then 
F z = T2V z 9 = sf sin (9)V z 9 and therefore the work done 
by the rotating DW-like field (or the SMF) on the spin- 
particle for one traversal (left to right) is 

ZR 7T 

W L ^ R = J F~ z dz = sfj sin (0)d0 = -2s f , (15) 



where ^l,r are the leftmost and the rightmost position 
of the spin-particle on the wire far from the region of 
spatial variation of the field. This result has been derived 
with the single assumption of adiabaticity. The adiabatic 
condition is s^ <C s and for it to hold it is necessary 
that the two components of averaged over the rapid 
precession, simultaneously satisfy the latter, i.e. 

s^l = 7\/jb = svV z 6/jb <C s , (16) 
s^ = %hb = sf cos (0)/7& < s . 

Thus, considering the maximum attainable values of the 
right-hand sides and since max(V 2 0) = l/z w , the neces- 
sary conditions for adiabaticity are 

V*w « /l, / « /l (17) 

where t w is the time it takes for the spin-particle to cross 
the DW-like region of width z w 9 and /l = 76 is its Lar- 
mor precession frequency about the field b. 

In order to mimic the typical strength of the exchange 
interactions (~ 1 eV) in our classical simulations, we have 
used 7 = 2/j,B/fi and b = 10 5 T. We have simulated an en- 
semble of TV = 700 noninter acting spin-particles, confined 
in a 400 A-long atomic wire. The particles start at ran- 
dom positions within two regions near both wire ends and 
with velocities identical in magnitude (vq = 8 A/fs) but 
random in sign. The DW-like region has z w = 5 A (see 9 ) 
which is similar to the z w , fitted to the relaxed DW pro- 
file in the quantum simulation (in atomic spacings) and 
the typical passage time is t w ~ 1 2 fs. We have used 
frequencies |/| < 0.2 fs -1 . Since our / L = 76 = 17.6 fs -1 , 



FIG. 6: (Color online) The spatial particle imbalance AN/N 
and the potential energy difference AW as a function of the 
angular frequency / of rotation of the DW-like field. The 
insets represent the change of the velocity and the longitudinal 
component of the spin of a particle for one left-to-right (red 
solid) and one right-to-left (blue dashed line) traversal of the 
wire for / = 0.2 fs -1 . 

these parameters well satisfy the adiabatic conditions of 
Eq. (17). 

We have integrated numerically Eqs. (11) and found a 
stationary difference in the number of particles to the left 
and to the right of the DW-like region, AN ca \ c = TVr — 
A/"l, developing in time and depending linearly on the 
frequency / of rotation of the field (see Fig. 6). By energy 
conservation, A7V ca i c converts to a potential energy shift 

Alcaic = 2m e v%AN calc /N (18) 

and the latter is a manifestation of the SMF work 
— Wl->r, derived in Eq. (15). Eq. (18), relating of par- 
ticle imbalance to SMF, is only valid if the particles have 
enough initial kinetic energy to traverse the wall from 
both sides, which, from one of the sides, means climbing 
the SMF ramp. Thus the requirement 

m e vl/2 > 2sf (19) 

sets a lower limit on the initial velocity of the spin- 
particles in our simulations, for a given /. 

Within the adiabatic regime the dependence of AW caL \ c 
on / is found to be linear with a slope of (0.643 ± 
0.012) eVfs (see the right-hand side scale of Fig. 6) 
and agrees with the analytical prediction of 2s = h = 
0.658 eVfs in Eq. (15). 

The directions of the SMF observed in the quantum 
and the classical simulations agree with the one set by 
Eq. (15), i.e. the SMF is opposite to the direction of 
the angular velocity of the DW rotation if the itinerant 
spins are aligned parallel to the local field. Note that 
with the choice of the band structure in our quantum 
simulations the effect is carried by the down-spin, so that 
for all other parameters being equal the sign of the SMF 
is opposite to that of the classical model. In general, the 
direction and magnitude of the voltage drop is found to 
scale with the Fermi level spin-polarization rj = (D^ — 
Di)/(D^ + D\), being -D-m) the spin-up(down) density 
of states at the Fermi level, and uo as Vl^r = —rjhcu/e, 



5 



where uj > corresponds to an anticlockwise rotation of 
the DW (spin) about the z-axis. For the half-metallic 
case studied here rj = — 1. 




FIG. 7: (Color online) Plots of the SMF as a function of 
the DW rotation frequency in units of the Larmor precession 
frequency of the itinerant spins about the local field. We 
present the case of dynamics away from the adiabatic regime 
for both the quantum [panel (a)] and the classical [panel (b)] 
simulations. Dashed green lines in both panels correspond to 
a slope of 1. 

A further similarity between the quantum and the 
classical simulations, pointing to the classical origin of 
the "quantum" SMF, is that the quantum effect relies 
strongly on the adiabatic conditions set by Eq. (17). As 
illustrated by Fig. 7(a) the effect dies out completely 
above the Larmor precession frequency uo\, = J/h of the 



exchange coupled spins for any choice of band filling. 
The threshold in the classical case below /l [see Fig. 
7(b)] is an artefact of the classical model and occurs at 
/ = f c = (m e VQ/2)/h as determined by Eq. (19). 



V. CONCLUSIONS 

In conclusion we have demonstrated computationally 
the presence of a spin-motive force in a quantum-classical 
system with a spatially and temporally dependent order 
parameter. Our SMF has the same magnitude as the one 
described by Barnes and Maekawa 1 . We have considered 
the more general case of an order parameter of mixed spin 
and orbital character in which case a measurable voltage 
drop across the system could indicate the presence of the 
SMF. We have also presented an analytical classical ar- 
gument for the mechanism of the SMF in the adiabatic 
regime. The latter is supported by purely classical simu- 
lations of particles with intrinsic angular momentum in a 
magnetic field with the spatial and temporal properties 
of the order parameter in the quantum case. The result 
is the same in magnitude SMF A</> = tiw/e, where the h 
factor comes from the magnitude s of the intrinsic angu- 
lar momenta, considered to represent the electron spin, 
i.e. 2s = h. 
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